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Abstract: We show in the example of a one-dimensional asymmetric exclusion process that 
stationary states of models with parallel dynamics may be written in a matrix product form. The 
corresponding algebra is quadratic and involves three different matrices. Using this formalism 
we prove previous conjectures for the equal-time correlation functions of the model. 
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During the last few years the study of one-dimensional reaction-diffusion models has been 
of increasing interest. These models describe dynamical processes far away from thermal equi- 
librium so that in general their stationary probability distribution cannot be derived from an 
energy function. Therefore different techniques are needed in order to determine the stationary 
properties. An exact method which turned out to be very successful is the so-called matrix prod- 
uct formalism Q-Q. This formalism can be seen as a generalization of stationary states with 
a product measure in which products of numbers are replaced by products of non-commutative 
algebraic objects. By representing these objects in terms of matrices, the stationary state and 
all equal time correlation functions can be derived exactly. Up to now this technique has been 
applied mostly to systems with sequential dynamics (continuous time evolution) where the sta- 
tionarity of the state is related to an additive cancelation mechanism from site to site. However, 
many systems, for example traffic models |§, are defined by parallel dynamics rather than se- 
quential updates. Therefore it is of interest to find applications of the matrix product technique 
to systems with parallel dynamics. The present work discusses this problem in the example of 
a one-dimensional asymmetric exclusion process with parallel updates ||. A modified matrix 
product formalism is presented in which a new multiplicative cancelation mechanism plays an 
essential role. The corresponding matrix algebra is derived and finite-dimensional representa- 
tions are given. This allows to prove exact results for the equal-time correlation functions which 
were already given as conjectures in Ref. Q. 

Let us first recall the matrix product formalism for one-dimensional reaction-diffusion models 
with sequential dynamics and open boundary conditions. A two-state model with L sites is said 
to have a matrix product ground state if the stationary probability distribution Po(t\, t%, . . . , t£) 
can be written as 

L 

Po(n,r 2 ,...,r L ) = Z- 1 (W\ HfrD + {l-r 3 )E)\V) (1) 

3=1 

where Tj € {0, 1} is the occupation number at site j. E and D are square matrices acting 
in an auxiliary space which may be either finite or infinite dimensional. The probabilities are 
the expectation values (W\ . . . \ V) of the matrix products normalized by the constant Z = 
(W\(D + E) L \V). Formally we may rewrite Eq. (|lj) tensor product 

|P ) = Z- l {W\ f^rV) (2) 



where the vector |Po) represents the stationary probability distribution in configuration space. 
The matrix product representation is a powerful tool since it allows various physical quantities 
like the particle density 

{W\C^DC L -i\V) (r -n + F\ ™ 

{t > )l - — WicW) — (3) 

to be computed directly. Higher correlation functions are given by similar expressions in which 
C plays the role of a transfer matrix. However, a special mechanism is needed in order to ensure 
that the state in Eq. (||) is indeed a stationary one. For models with sequential dynamics this 
mechanism amounts in an additive cancelation from site to site: Assume that the time evolution 
of the system is described by a master equation 4i\P) = —H\P) with a time evolution operator 

H = J2j=i + + ' w here is a 4 x 4 interaction matrix and and are 

2x2 matrices for particle input and output at the ends of the chain. Then the matrices E and 
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D have to be chosen such that the application of the interaction matrix hjj+i to a pair of sites 
results in a local divergence-like term on the right hand side 



E\ „ E 



DJ®\D 



b) » (I) " (d) 8 (6 



(4) 



where i? and D are again matrices in the auxiliary space. By summing up the two-particle 
interactions, all these contributions cancel in the bulk of the chain. The remaining terms at the 
boundaries have to be canceled by the surface fields for particle input and output: 

< W\ h m (£)=_< wl (?) , „<«> (g) IV >- (!) |V > . (5) 

The simplest case E = D = and its generalization to spin one problems has been considered 
in Ref. ||. Another system which has been investigated in detail is the (asymmetric) exclusion 
process where E = -D = 1 g. In both cases one is led to an quadratic algebra of two objects E 
and D Q. In a similar way matrix product ground states were found for particular three-state 
models [||. Also excited states can be described with a matrix ansatz || where E = —D has to 
be chosen as a time-dependent matrix leading to a quadratic algebra of three different objects. 
Taking E and D as independent matrices, it was also possible to find the stationary state of 
particular models with particle reactions [Q]. 

So far, the interest has been focused mainly on stochastic models with continuous time 
evolution. However, similar techniques can be used for systems with parallel dynamics. A first 
example of this type was given in Ref. where the transfer matrix for a deterministic model 
of directed animals on a strip was investigated. It is the aim of the present work to point out 
that there could be a broad spectrum of applications to reaction diffusion models with parallel 
dynamics. For this purpose we consider a one-dimensional asymmetric exclusion process with 
parallel updates which was originally introduced by G. M. Schiitz in Ref. ||. In this model 
particles move on a one-dimensional lattice with L = 2N sites and open boundaries. The bulk 
dynamics is deterministic and consists of two half time steps. In the first half time step particles 
at odd positions move one step to the right provided that the neighboured site to the right is 
empty. Then in the second half time step the particles at even positions move to the right in 
the same way. In addition particles are injected (removed) stochastically with rate a (/?) at the 
left (right) boundary: 



a \ r\ c\ c\ a 

first half time step: o o o o o o o— — o P 

1 2N 



r\ 

second half time step: oooo««« o 

1 2N 



The corresponding transfer matrix therefore consists of two factors T = T^Fx 

Ti = C®T ® ...®T ®K = £®T® {N -^ ®K (6) 
T 2 = T ®T ® ...®T = T® N 
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where T, C and 1Z are the matrices for hopping, particle input and output: 



T 



( 1 \ 

110 



\ 1 J 




n 




(7) 



The phase diagram of this model shows two phases. For a < (3 the system is in a low-density 
phase with an average particle density p = a/2 < 1/2 whereas in the high-density phase a > (3 
one has p = 1 — (3/2 > 1/2. The total current in the thermodynamic limit is given by j = 
min(a, (3). The physical behaviour is closely related to that of asymmetric exclusion models 
with continuous time evolution |J (there is an additional phase with maximal density in the 
latter case). It plays a role in traffic models |8| as well as in polymer physics [10]. Related 



models with deterministic dynamics can be found in Ref. [11] and the influence of defects has 
been studied in Ref. JL^]. 

As we are going to show below, the stationary state of the exclusion model (|6j)-([5) c & n be 
written as a matrix product state with alternating pairs of matrices (E, D) and (E, D) so that 
the probability to find the system in the configuration (ri, t 2 . . . t 2 n) is given by 

N 

Pq(ti,t 2 ,...,t 2N ) = Z- X (W\J[ [{r 2i -ib + (l-r 2i -i)£) (r 2i D + (l-r 2i )^)] \V) . (8) 

i=i 



As in Eq. (|2|), we may also write 

|P ) = Z- 1 {W\{f ) )^ KD 



(!) 



E D )W) 



(9) 



z-\w\ 



(I) 



\V) 



where Z = (W\ \^(E + D) ®(E + D)J \V). It is obvious that in this case the basic mechanism 
leading to a stationary state has to be different from the usual one for continuous time evolution 
operators. Instead of the additive cancelation from site to site we now need a multiplicative 
mechanism suitable for stationary states states of parallel transfer matrices T\Pq) = \Pq)- In the 
case of the above exclusion model this mechanism turns out to be very simple. Let us assume 
that in each time step the two pairs of matrices (E, D) and (E, D) are exchanged: 



Ti\Po) 



(W\ 



(d) - (!) 



E- 



\V) 



T\P ) 



T 2 Ti|P ) 



This exchange mechanism can be realized by 

'E\ „ (E^ 



T 



(10) 
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which is equivalent to the algebra 



n( E D )\v) = (%)\v) 



E- 



[E,E] 
ED 
DE 



[D,D] 
[E,D] 




(12) 
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and the boundary conditions 



(W\E(l-a) = {W\E (1-P)D\V) = D\V) 

(W\(aE + D) = (W\D (E + f3D)\V) = E\V) 

The commutation relations (|l^) involve four different matrices. However, only three of them 
are independent since the matrix product in Eq. (||) is invariant under the transformation 

E^U~ l E, D^U~ l D, E->EU, D -» DU . (14) 

Because of [E+D, E + D] = it is possible to choose a basis in which both operators E + D and 
E + D are diagonal. Taking now U = (E + D) 1 / 2 ^ + D)~ l l 2 both operators become identical 
so that we may add the relation 

C = E + D = E + D. (15) 

Eliminating E and E we therefore obtain a quadratic algebra of three independent objects 
defined by three bulk equations 

DC = DD = DD, [D — D,C] = (16) 

and two boundary relations 

(W\ (D - aC - (1 - a)£») = , ((1 - p)D - t)) \V) = . (17) 

Matrix product states based on quadratic algebras with three objects were first studied in 
Ref. |§. A detailed analysis of algebras with more than two objects and their representations 
will be given in Ref. [13|. 

In order to check the consistency of the algebra (|l~5|)-(16) let us show that the expectation 
value of any sequence of operators is given uniquely on an abstract level. For only two operators 
it can be verified by hand that 



(W0CV) _ (W\DD\V) _ W -^Z£L- W{W \CC\V) (18) 



w^iv) = {a ^ a m-v+P mv) <19> 

so that any expectation value of two operators is a given number times Z = (W\CC\V) . In order 
to check the consistency of the algebra for products of arbitrary length, it is more convenient to 
use a different basis of operators which is defined by the invertible transformation 



X = —(D-aC+(a-l)D 



Y = ±-((l-p)D-D) (20) 



a(5 



In this basis, the bulk algebra ( p^q ) reads 

[X,S] = [Y,S] = 0, YX = {l-a)SY + (l-(3)XS-(l-a)(l-(3)S 2 (21) 
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and the boundary relations (|l7|) become particularly simple: 



(W\X = 0, Y\V) = 0. (22) 



As can be seen easily, the application of the bulk relations ( |2l[) allows every product of 2N 
matrices X, Y and S to be ordered as a linear combination of terms like x n S 2N ~ n ~ m Y m . 
Since the only nonzero expectation values of these terms is (H /r |S ,2Ar |y), the expectation value 
of any product of 2A^ matrices is a well-defined number times (WIS^IV). The actual value 
of (VF|(S' 2Ar |y) is irrelevant since it is canceled by the normalization constant Z = (W\C 2N \V). 



Thus it is obvious that the algebra (21)-(|22|) determines the ground state |Po) uniquely on an 
abstract level. We should emphasize that the mathematical structure of this algebra is different 
from that for exclusion models with continuous time evolution where one has linear terms in 
the bulk algebra (e.g. DE = D + E). Whereas in the latter case any expectation value can be 
reduced to the empty bracket (W|V), the algebra ( |l6|) does not allow to reduce the number of 
factors in a matrix product. Instead of this we have shown that by means of the algebraic rules 
the expectation values of all words with the same number of factors are linearly dependent. 



The algebra (|lq)-(|17j) can be represented by two-dimensional matrices. For a ^ f3 a repre- 
sentation in which C is diagonal is given by 

°) D 1= ( » ° ) I) (23) 



-a/3 

(Wi\ = (a, 1-a), |Vi) = rlf^j (24) 

The normalization constant in this representation can be computed easily and reads 

Zi = (1 - (3) a 2N+l - (1-a) (3 2N+1 . (25) 

As already mentioned, the matrix C acts like a transfer matrix between the points of the 
correlation functions. Therefore the length scales to be expected are essentially given by the 
quotients of the eigenvalues of C . Thus in the present case the correlation functions involve 
only a single length scale, namely log(a//3). This length scale diverges at the phase transition 
line a = f3 where the constant Z\ vanishes so that the above representation becomes singular. 
It turns out that in this case the operator C cannot be diagonalized so that one has to use a 
different representation where C has a Jordan normal form: 




r '-'-f " 1 ) D2 ={0 1 ) D2 ={0 1-a) ^ 

(W 2 \ = (1, 0), \V 2 ) = (/„.) (27) 



Because of 

<* - ( J J ) < 28 > 

the normalization constant Z is now linear in the system size: 

Z 2 = 1 + 2N(1 - a) . (29) 
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Using the matrix product formalism it is now easy to derive explicit expressions for the 
equal-time correlation functions. Following the ideas of Ref. Q, we first compute the n-point 
functions of the operators 

_ T2j - a _ T 2 j-1 , . 

m 3 - mj-i - ^ (30) 

Denoting the corresponding matrices by 

j (D - aC)/ (1 - a) if j even 
] \ D/(l - (3) if j odd 1 j 

and assuming that the positions ji ■ ■ ■ j n are are chosen in increasing order these correlation 
functions are given by 

(VhVh ■ ■ ■ V jn ) = \ {W\C^~ l F n C**-*- 1 F n CI- ' 1 . . . C^-i- 1 F jn C 2N ~ jn \V) . (32) 



Using the representations (p3])-(|29|) it is easy to check that 

Fj C k -*- x F k = FjC k -i (j < k) (33) 

so that the n-point correlation functions reduce to the one-point function (rjj): 

= (Vn) = Z- x {W\C^ l F n C 2N -^\V). (34) 
For a ^ (5 the one-point function reads 

( V2j ) = ±-a™+^(l-P)(a^-^) (35) 

faw-i> = ±oF+W(aF-Hl-P)-lP-\l-a)) (36) 
whereas at the transition line a = (3 we have 

(V2j) = ^2j(l-a) (37) 

(V2j-i) = ^(a + (2j-l)(l-a)). (38) 



Z 2 

Although we used the two-dimensional matrices at this point, Eqs. (|34[)-(|37|) do not depend on 
the choice of the representation since we have shown that the expectation value of any sequence 
of operators is uniquely given by the commutation relations of the algebra. 

Resubstituting Tj into Eq. (34) we obtain an exact expression for the n-point density corre- 
lation functions (t j 1 t j 2 ■ ■ - Tj n ). Denoting Oj = j mod 2, they read 

n n k—l 

(j h r j2 ...r jn ) = a n Y[ai + £ (j[ 07) ^(l + aa k - - /3<r fe ) (r, jk ) . (39) 

i=l k=l i=l 

As a special case this formula includes the two-point correlation functions which have been given 
as conjectures in Ref. |J. 

The example of the asymmetric exclusion model shows that the powerful matrix product 
formalism can be applied successfully to models with parallel dynamics. Since models of this 
type are widely studied, it would be interesting to find further examples in order to understand 
under which conditions the matrix product technique can be applied. In particular it would 
be interesting solve the same model on a ring in the presence of a defect. From this one could 
learn how to solve the full exclusion process (with stochastic hopping in both directions) on a 
ring with a defect [ 14 ] . Despite of intensive efforts the exact solution to this problem is not yet 
known. 
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